Avoided level crossing statistics in open chaotic billiards 
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We investigate a two-level model with a large number of open decay channels in order to describe avoided 
level crossing statistics in open chaotic billiards. This model allows us to describe the fundamental changes of 
the probability distribution of the avoided level crossings compared with the closed case. Explicit expressions 
are derived for systems with preserved and broken Time Reversal Symmetry (TRS). We find that the decay 
process induces a modification at small spacings of the probability distribution of the avoided level crossings 
due to an attraction of the resonances. The theoretical predictions are in complete agreement with the recent 
experimental results of Dietz et al. (Phys. Rev. E 73 (2006) 035201). 
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It is by now established, that classical chaos manifests it- 
self in universal spectral fluctuation properties of the eigen- 
values of the corresponding quantum system. They coincide 
with those of random matrices from the Gaussian Orthogo- 
nal Ensemble (GOE), if Time Reversal Symmetry TRS holds, 
from the Gaussian Unitary Ensemble (GUE), if TRS is broken 
fliii. 

Investigations of the universality of spectral fluctuation 
properties of classically chaotic systems range from nuclear 
physics fio, to systems in other areas, like microwave bil- 
liards fTllsiMfloll. optical experimen ts llTl[[l2ll . quantum dots 
JIUHL and acoustic setups lfl5l[l6lll7ll . In systems depend- 
ing on a global parameter, the correlations between eigenval- 
ues at different parameter values show a universal behavior, 
which a gain is well described by random matrix theory (RMT) 
ifisi [T9II20I I2TI1 . In some cases, implying a local parameter, 
RMT fails, as reported in reference ll22[l . 

In an experiment presented in ll23ll the spectral properties of 
a superconducting microwave billiard, whose boundary was 
varied parametrically, were investigated. It models a quantum 
billiard of corresponding shape, whose classical dynamics is 
chaotic. The observed deviations from the expected GOE be- 
havior were attributed to the measurement process. Indeed, 
resonance spectra of a microwave billiard are measured by 
connecting it to the exterior via emitting and receiving anten- 
nas. Thus the resonator is an open system with the antennas 
acting as single scattering channels. The influence of the flux 
of microwave power flowing from the emitting to the receiv- 
ing antenna on the spectral properties of the system is so weak 
that it cannot be detected through spectral measures like the 
nearest neighbor spacing distribution or the statistics at a 
fixed value of the parameter The distribution of the avoided 
crossings of the eigenvalues as function of the parameter on 
the other hand showed deviations from the GOE result, which 
were attributed to the openness of the resonator. These as- 
sumptions were confirmed by numerical simulations based on 
a random matrix model for parameter dependent, chaotic and 
open systems. The aim of the present paper is the derivation of 
an analytic expression for the avoided-crossings distribution 
of such systems. It goes in line with that for the corresponding 



distribution for closed systems with and without TRS, which 
is based on an ensemble of two-dimensional random matrices 
M- 

We develop our approach within the framework of an ef- 
fective Hamlitonian model [[24|l. To describe statistical prop- 
erties of avoided crossings in open systems, we introduce the 
effective Hamiltonian i?e// which depends on a continuous 
parameter /i through its Hermitian part ll23ll 

H,ff{fi)^H{^i)-^VV^ . (I) 

Here H{fi) is the Hamiltonian of the closed system modeled 
by a 2 X 2 random matrix and iVV'^ /2 is an imaginary po- 
tential describing the coupling to the environment in terms of 
M open channels. The 2 x AI matrix V contains the coupling 
amplitudes y,™ which couple the nth level to the mth open 
channel. As a result, the eigenvalues of the effective Hamilto- 
nian are complex, e± = E± — ^T±, where Ej- and T± are, re- 
spectively, the two eigenenergies and the two spectral widths 
of the 2-level model. For the study of statistical properties 
H is replaced by a Gaussian random matrix Izsll and the ma- 
trix elements y,™ are chosen to be Gaussian-distributed with 
zero mean and variance = 2AA, where A is the coupling 
strength and A is the mean level spacing of the closed system 
1E6I1 . In the eigenbasis of H{jj.) (the ^-dependence is omitted 
in the following), the effective Hamiltonian is written as 
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where £^12 are the //-dependent eigenenergies of H {E2 > Ei 

is assumed) and r„p = X]c=i KI"^™- Note that the model is 
applicable only as long as the coupling is weak enough so that 
the spectral widths remain of the same order of magnitude ll27l 
I28I1 . The complex eigenvalues of the effective Hamiltonian 
Hef f in Eq. ^ read 

Ei+E2-^iTn+T22)±VD 
£± = n 
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with 

D = - E^) + ^{T22 - Tii))^ - riaFai . (4) 

The spacing between the two eigenergies, d — e+ — e- can be 
read off from (O and 



d = Re(/D) . 



(5) 



Considering the limit of a large number of open channels M 
in the weak coupling regime, we may apply the central limit 
theorem and replace the random variables depending on the 
coupling amplitudes by their averages. 
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Then the spacing d is given by 

Vs2 -Mcr4 ifs>VMcr2 
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where s ~ E2 — Ei is the spacing of the eigenenergies of 
the closed system. Note that Ma" = var(r)/2 implies that 
the modifications on the spacings due to the openness of the 
system are related to the fluctuations of the spectral widths 
Q. In the hmit M ^ 00 and cr^ ^ with Ma"^ = (T) 
fixed, var(r) and thus the spacing between eigenenergies 
of the open system converges to that of the closed system, 
d ^ s, in spite of non-vanishing losses. One of the effects of 
the imaginary potential is that the eigenvalues mutually attract 
each other along the real axis fsoll. As this attraction increases 
when s decreases, the local minima of both d and s coincide. 
In other words, the values of the parameter ^ at the avoided 
crossings are the same for the closed and the open system. 
Accordingly, in the derivation of the distribution of avoided 
crossings c of the open system, the spacings s are assumed to 
be distributed as the avoided crossings of the corresponding 
closed system. With Eq. (ITJ the probability distribution of the 
avoided level crossings p{c) is given by 

p{c) = (5{c)e{\fMa^ - s)^ 

+ (5{c- ^/ s"^ - Ma"^e{s - ^la^)^ , (8) 

where 9 is the Heaviside step function and the triangular 
brackets denote averaging with respect to the spacing s. 

For closed chaotic systems with TRS the probability dis- 
tribution of avoided crossings has been calculated by Za- 
krzewski and Kus ifisll . 
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where the mean value of s is given by (s) = a^JYfn. Aver- 
aging over s yields 



p{c) = erf 



Ma' 



5{c) 
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where a fixes the average of c. Note that the behavior of 
p{c) at small spacings differs strongly from the GOE predic- 
tion (|9]l. The linear behavior of the distribution induces a dip 
and the local minima of the spacings have a zero-crossings 
contribution leading to the presence of a (5-peak at the origin. 
This peak is neither restricted to 2-level models nor to a large 
number of channels. It was also found numerically in ll23ll 
where an effective Hamiltonian with 1000 levels and M = 3 
open channels was considered and is characteristic of non- 
Hermitian Hamiltonians of the form Eq. © IHIllllllll. 

In the experimental setup ll23ll . only finite spacings could be 
measured due to the discrete sampling of the data. Therefore, 
to compare theory and experiment it is more convenient to 
consider the distribution of non zero avoided crossings c'. 



p{c') 
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The analysis can be extended to open, parameter depen- 
dent, chaotic systems with broken TRS. For closed systems of 
this type the probability distribution of avoided crossings ifisll 
reads 



p(s) 



(12) 



Using the same effective Hamiltonian model Q i.e. consid- 
ering real coupling amplitudes, the probability distribution of 
avoided level crossings is derived using (|8]l and ( fT2] i. 



- Cl - p-M^V{4«') 



V{c) = (1 - e 



)<5(c) 



(13) 

Again a 5-peak appears at c = due to the attraction of the 
eigenvalues on the real axis. However, in contrast to the GOE 
case, the probability distribution of the non zero avoided level 
crossings coincides with that of the closed system given in 
Eq. (HI 



p{c!) 



2a2 
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This robustness of GUE was previously observed in room 
temperature microwave billiards with broken TRS fl35l1 . 

To analyze the evolution of both distributions for a small 
or not too large number of channels, numerical random ma- 
trix simulations were performed, where the eigenvalues of the 
closed, parameter dependent system were chosen as the eigen- 
values of the random matrix 



_ff (/i) = R\ cos /i + ^2 sin /i , 



(15) 



simulating the closed system (see II23I1 ). Here H\ and Hi be- 
long to the GOE or the GUE for the simulation of systems 
with or without TRS, respectively, and the coupling ampli- 
tudes with n = 1, • • • iN and m = 1, • • • , M are ran- 
dom Gaussian variables. Note that this model ensures that 
the mean level spacing is independent of /i ll20ll . In the sim- 
ulations, the matrices El\ and are of size 1000 x 1000, 
the variances of their elements are chosen equal and such that 
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A = 1/1000, the parameter /i G [0, 7r[ is discretized in steps 
of Sfi = 7r/1300. To ensure a fairly constant mean level 
spacing only the 400 resonances at the center of the Wigner 
semicircle were kept. To mimic the experimental resolution, 
a cut-off Co = O.lA is introduced such that only values of c 
larger than cq are used to build the numerical distributions. 
Furthermore the average of the spectral widths is fixed to 
(r) = 0.5A, well away from the strong coupling regime ll23ll . 
A comparison between analytic and numerical results is pre- 
sented in FiglT] 



GOE P(c') GUE 
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FIG. 1: (color online) Probability distributions of non zero avoided 
level crossings for GOE (left side) and for GUE (right side). The 
number of open channels is M=l, 3, 5, 10, 20 and the coupling 
strength A = a^/{2A) equals A = 0.250, 0.083, 0.05, 0.025, 0.013 
from top to bottom. The histograms show the numerical simulations, 
the analytic distributions are shown as straight lines. The dashed 
curves result from a fit of Eq. dllb with A as a parameter to the 
numerical distributions, resulting in the effective coupling strengths 
Xeff = 0.020, 0.038, 0.036, 0.028, 0.016. For all curves, a is cho- 
sen such that (c') — 1. 



For the GOE case a good agreement between the numer- 
ical and analytical descriptions is found for A/ > 5. For a 
smaller number of channels, the histograms are reproduced 
by choosing the coupling strength A = cr^/(2A) in Eq. ( fTTI) 
as a parameter to obtain an effective coupling strength Ae/ / 
by means of a fit based on a least square algorithm. Thus, 
it appears that the expression ( fTTI ). derived using the central 
limit theorem, can be extended to any M, considering A as a 
free parameter. 

The right column of Fig. [1] shows the results obtained for 
the GUE case. The prediction is in excellent agreement with 
the numerical results except for the case M = 1. This is 
due to the small number of events at small distances of p{s). 
Indeed, while important changes appear for the GOE due to 
the large number of small avoided crossings, the GUE case 
is only slightly modified because of a vanishing density of 
avoided crossings at the origin for the closed case. Note that 
the distribution ( fT4l i is independent of the coupling strength 
such that a fitting procedure is not possible. 




FIG. 2: (color online) In boxes the experimental distribution of 
avoided crossings I23I1 . The continuous line shows the analytical 
prediction with A=0.058 (obtained through a least square procedure). 
The vertical bars represent the numerically obtained distribution with 
M = 3 and A=0.02 123il . For all curves, the average is chosen such 
that (c') = 1. 

Now, let us finally compare the analytical prediction (fTTI ) 
with the experimental results of ll23ll obtained using a su- 
perconducting microwave cavity, thus minimizing dissipative 
processes. Three antennas were attached to the cavity: they 
correspond, in our model, to M = 3 open channels ll36ll . 
Absorption into the walls could be mimicked by additional 
fictitious weakly coupled channels ifTol l26l [stIi . however, its 
influence can be safely neglected in the analysis of the ex- 
perimental data. Due to the lack of an analytic expression 
for the distribution of avoided crossings in open systems, the 
experimental distribution was compared with numerical sim- 
ulations based on an effective Hamiltonian of the form Eq. ([T]) 
with the parameter dependent Hamiltonian of the closed sys- 
tem given in Eq. ( fTSt . In reference ll23ll . a good agreement 



4 



beween both distributions was shown for values of A ~ 0.02. 
Note that due to the lack of an analytical expression this value 
was not determined from a fit and thus is only vague. A com- 
parison between the experimental (boxes) and numerical his- 
tograms (vertical bars) is shown in Fig.|2l The numerical dis- 
tribution has been computed with A = 0.02. The analytic 
result obtained through a fit using A as a parameter (contin- 
uous curve) is also shown in Fig. |2] The 2-level model re- 
sult follows closely the experimental histogram. This con- 
firms the interpretation drawn in ll23ll that the deviation of the 
avoided-crossings distribution from the predicted GOE result 
for closed systems is due to the measurement process, i.e. the 
influence of the three antennas, which couple the resonator 
modes inside the resonator to the exterior. 

In summary, we have derived an analytic expression for 
the distribution of the avoided crossings of the resonances of 
quantum chaotic open systems based on a simple two-level 
random matrix model. Analytical results prove that the open- 
ness essentially modifies the avoided crossing distributions at 
small spacings. The theoretical predictions are in excellent 
agreement with numerical random matrix simulations for a 



number of open channel M > 5 in the GOE case and for 
A/ > 3 for GUE systems. For systems preserving TRS with 
a small number of open channels good agreement is achieved 
by using the coupling strength A as a fit parameter. 

Finally, let us mention that the 2-level model can also be 
used to calculate the nearest level spacing distribution (NLSD) 
for open chaotic systems with a large number of channels. 
Whereas, for the strong coupling regime, the NLSD is sub- 
stantially modified by the introduction of dissipation Issll . in 
the weak coupling case, due to a vanishing density at small 
spacings for both GOE and GUE, the NLSD for open systems 
will only be moderately modified, as pointed out above. 
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